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ON AN ELLIPTIC EQUATION ARISING EROM COMPOSITE 

MATERIALS 

HONGJIE DONG AND HONG ZHANG 


Abstract. In this paper, we derive an interior Schauder estimate for the 
divergence form elliptic equation 

Di{a{x)Diu) = Difi 

in where a{x) and fi{x) are piecewise Holder continuous in a domain 
containing two touching balls as subdomains. When fi = 0 and a is piecewise 
constant, we prove that u is piecewise smooth with bounded derivatives. This 
completely answers a question raised by Li and Vogelius in dimension 2. 


1. Introduction 

Ill this article, we consider second-order divergence type elliptic equations with 
discontinuous coefficients and data 

,r2^ ■— — ^ifi 1 ^ i (^'^) 

where P is a bounded subset of 

a{x) = aoXBri(0,ri)UB.,2(0,-r2) + XR2\(B..j(0,ri)UB^2(0.-r2))) 
ao > 0 is a constant, ri, r2 € ( 0 , oo), and x is the indicator function. This problem 
was raised by Bonnetier and Vogelius [ 5 ], and can be considered as a simplified 
model for composite media with closely spaced interfacial boundaries. Here T) 
models the cross-section of a fiber-reinforced composite and the balls 5 ^( 0 , ri) 
and Br2(0, —r2) represent the cross-sections of the fibers; the remaining subdomain 
represents the matrix surrounding the fibers. Moreover, a{x) is the shear modulus, 
which is a constant on the fibers, and a different constant on the matrix surrounding 
the fibers. The function u stands for the out of plane elastic displacement. 

Elliptic equations and systems arising from elasticity have been studied by many 
authors. See, for instance, [iEiiiiisiiniiTiiiaiii]- In 0 , Chipot, Kinderlehrer, 
and Vergara-Caffarelli considered divergence type uniformly elliptic systems in a 
domain D C consisting of finite numbers of linearly elastic, homogeneous, par¬ 
allel laminae, which models the equilibrium problem of linear laminates. In [ 9 ], 
Li and Vogelius studied divergence type elliptic equations in a bounded domain 
V C R'^, where V can be divided into finite numbers of subdomains with (7^’“ 
boundaries. The coefficients of the equations and data are Holder continuous in 
each subdomain up to the boundary, but may have jump discontinuities across the 
boundaries of the subdomains. Under these conditions, they proved a global W^’°° 
estimate and a piecewise estimate, for any /3 < Notably, their estimate 

does not depend on the distance of discontinuous surfaces, which indicates that by 
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Figure 1. 

an approximation argument, interfaces may touch each other, e.g., the geometry 
shown in Figure[T] Later, Li and Nirenberg [TU] extended the result in to elliptic 
systems under the same condition. They were able to improve the piecewise 
estimate in to any fj € (0, 

Regarding the operator in (HD, Bonnetier and Vogelius [5] first considered the 
Dirichlet value boundary with ri = r 2 = 1 and fi = 0: Lipu = 0 in (D and u = (j) on 
dV. They showed a global regularity result that the solution u G Later, 

Li and Vogelius [9] extended the result in [5] and proved that when ri = 7-2 = 1, 
fi = 0, and V = with i?o sufficiently large, the weak solution u is piecewise 
smooth, i.e., 

nG C7°°(Ri(0,l)), uG C°°(Ri(0,-l)), n G \ (Ri(0,1) U Bi(0,-1))), 

where K is any compact subset of Bfi^. Then they asked the following natural 
question: can we drop the condition that i?o being sufficiently large? 

In our first result, we answer this question by proving that i?o > 2 is sufficient 
to guarantee that u is piecewise smooth in the interior of Br^ . 

Theorem 1.1. Let Rq > 2 and g G H^^^{dBRg). Suppose u is a weak solution of 

Di{a{x)Diu) = 0 in Br^, u = g on BBr^, 

where 

a{x) = aoXBi(o.i)uBi(o,-i) + Xbb(,\(Bi(o.i)uBi(o,-i))- (1-2) 

Then 

uG C°“(Ri(0,l)), uG C°“(Ri(0,-l)), uGC°“(X\(Ri(0,1)URi(0,-1))) 

for any compact set K C Br^ . 

To prove Theorem 11.11 we borrow some ideas from [9]. In [9], Li and Vogelius 
constructed a sequence of piecewise smooth solutions {uj} to 

Di{a{x)Diu) =0 in 

the linear combinations of which are dense in H^y^^dBRg) for Rq sufficiently large, 
where H^y^{dBRg) denotes the space of functions even in xi with finite norm 
for s > 0. The precise definition can be found in Section 2. Therefore, the solution 
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Figure 2. 

u to the Dirichlet problem with the boundary condition u = (j) £ 
can be approximated by linear combinations of m^’s. Hence, by a classical elliptic 
regularity argument, one can show that \D^u\ < oo in each subdomain for any 
fc > 0. 

In this paper, we carry out a more careful analysis on {uj} to show that i?o > 2 is 
sufficient to guarantee that {uj} forms a Schauder basis for Precisely, 

it is obvious that 

{ej,j > 0} := {(-1)-^ cos(2j6»), sin((2j + 1)0) J > O} 

is an orthogonal basis of H^y^^dBn^). Each Uj can be written as a linear combina¬ 
tion of Cj’s, i.e., Uj = Mj^k^k- We show that the infinite dimensional matrix 

M := {Mkj)^j^Q define a bounded and invertible operator on a Hilbert space P. 
For the definition of I®, see Section 2. An important observation in our proof is 
that the submatrix is diagonally dominant by column. From this, we 

deduce that the map induced by M is invertible, which further implies that {uj}j>o 
forms a Schauder basis of = i?o}). The remaining proof then follows the 

lines in [3]. 

Another natural question to ask is if the geometry of the domain where the 
equation is satisfied affects the smoothness of the solution around the origin? In 
other words, if Lu = 0 in T>, is it necessary that 1) contains a ball with radius Rq > 2 
for u to be piecewise smooth around the origin? Or T? can be any neighborhood of 
the origin? Our second result answers this question by proving an interior Schauder 
estimate for the non-homogeneous equation cm in a general domain. Furthermore, 
we break the symmetry of the coefficients, meaning that a{x) can be two different 
positive constants ag and bg in the two balls with different radii. 

Theorem 1.2. Let ri,r 2 G (0, oo), 7 G (0,1), and n > 0 be an integer. Assume 
that B (ZM? is a bounded open set. Suppose that for any i, fi is piecewise C"’'*', 
i.e., 

/, e n Br, (0,n)), /, e n H,,(0, -r2)), 

and 

h G (H,,(0,ri)UH,,(0,-r2))). 
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Let u be a weak solution to 

D^{a{x)Diu{x)) = Dif^ inV, 

where 

a{x) = aoXSrJO.r-i) + boXBr^^O-r^i) + XR 2 \(Bri( 0 ,r'i)UBr 2 ( 0 ,-i- 2 )) 

Define = {x £ D, dist(x^ dD) > e} for any e > 0. Then 

u£C^+^’^V,nBrAO,ri)), u£C^+^’^V,nBr,{0,-r2)), 

and 

u £ \ (B,,(0,ri) U B,M-r 2 ))). 

In particular, when fi = 0, u is piecewise smooth in up to the boundary. 

For the proof, first we find a conformal mapping which maps two balls with 
different radii to two balls with the same radius, so it is sufficient to consider 
ri = r 2 and we denote the elliptic operator in (1131) with ri = r 2 = 1 by L. 
Then the conformal mapping F : z —>■ i/z maps {\z — \y\ = 1} and {\z + iy\ = 1} 
to {Rez = i} and {Rez = — respectively, where i = \/^ is the imaginary 
unit. We are able to construct Green’s function G{x,y) of the elliptic operator 
Lu = Di{A{x)Diu), where 

A{x) = aoX{^^> 11 + boX{^,<-L} + X{|xi|<i}- 

With the help of F and G{x,y), we obtain Green’s function G{x,y) of the elliptic 
operator L in which can be written as an infinite series of logarithmic function 
composed with smooth functions, for example, when x £ R^\ (i?i(0, l)Ui?i(0, —1)), 
and y £ Bi(0, 1), 

OO OO 

G(x, y) = Cl log \X_ 2 k{x) - y\ - C 2 ^(a/3)''"^ log \X 2 k-i{x) - y\, 

k=0 k=l 

where ci,C 2 ,a, and f are constants with |a|, \j3\ < l,y = fyi, — 2 / 2 ), and {Xk} are 
conformal maps and Xo{x) = x. Note that log|x — y\ is Green’s function of the 
Laplacian up to a factor. This observation allows us to implement some known 
results of the Laplace equation with piecewise Holder continuous data on the right- 
hand side. More precisely, the original problem is decomposed to understand the 
regularity of solutions to the following equations 

Am = A(/iXBi(o,i)), 

Au = A(/iXBi(o- 1 )), 

Am = A(/iXR2\(Bi(o,i)uBi(o,-i))), 

where in each subdomain fi £ C'”’''". By locally flattening the boundary, the first 
two equations can be further reduced to the case that fi £ in two half spaces, 
i.e., {x 2 > 0} and {x 2 < 0}. The detail can be found in the proof of Theorem 14.91 
Gase 1. The last equation needs an extension result to be reduced to the previous 
case. See Lemma [2.11 Gombining with the smoothness of each Xk, we are able to 
estimate all the derivatives of the solution. 

By a standard perturbation argument, we have the following corollary. 
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Corollary 1.3. Let 0 < A < A < c» and n > 0 be an integer. Assume a{x) is 
piecewise C"’'*', i.e., 

a{x) e C"’^(B,,(0,ri)), a(i) e (0,-r 2 )), 

and 

a{x) G (B,,(0,ri)UB,,(0,-r2))), 

and satisfies the ellipticity condition X < a < A. Suppose that for each i, fi is 
piecewise C"’'*', i.e., 

u G n Br, (0,n)), /, G n s,,(o, -r 2 )), 

and 

u G (i?,,(0,ri)U 5^,(0,-r 2 ))). 

Let u be a weak solution to 

Di{a{x)Diu) = D^fi 

in "D G Denote He = {a; G H, dist(x,dD) > e} for e > 0. Then u is piecewise 
(jn+i,'r jyj 2?e up to the boundary. 

This paper is organized as follows. In the next section, we introduce some nota¬ 
tion and preliminary results, which are needed in the proof of our main theorems. 
In Section 3, we prove Theorem ll.il In Section 4, we make necessary preparations 
and prove Theorem 11.21 and Corollary 11.31 


2. Notation and preliminary results 


In this section, we first introduce some notation used throughout this paper. The 
Einstein summation convention is applied in this paper. We use Bji{x) to denote 
the Euclidean ball in with radius R and center x. For simplicity, Bi{0, 1) and 
Bi{0, —1) are denoted by *Bi and ® 2 , respectively, and \ U $2 = ®o- When 
there is no confusion, we use Bii„ to denote the ball with radius Rq and center 
(0, 0). We use L to denote the operator Lri^r 2 when = r 2 = 1. 

Let H be a subset of and /3 G (0,1]. For any function /, we define 


[/]/3;X) = sup 


\f{x) - f{y)\ 

|a; - y\P 


: x,y G'D,x y 


and 

ll/ll/3;T> = II/I|Loo;T> + [/]/3;X>- 

We denote the space corresponding to || • Wp-xi by C^{'D). For nonnegative integer 
m, we dehne 

||/||m./3;T> = |1/|1 l„,;T> + [-D™/]/3;X>. 

The space corresponding to |1 • ||m,/ 3 ;i 5 IS denoted by 

Denote Ll^^(^{x : |a;| = i?o}) to be the set of real-valued functions on the 
circle {|a;| = Rq} which are even with respect to a:i. We use a similar notation for 
the Sobolev spaces Hgy^[{x : |a;| = Rq}) for s > 0. Note that L‘^y^[{x : |a:| = 
A-o}) = : 1x1 = i?o}). 

Let E be a Banach space over R. We say that a sequence {&„} in E is a Schauder 
basis of V if for every v € V, there exists a unique sequence {a„} of scalars such 
that 

00 

V — ^ ^ ^TL^ni 
n—{) 




6 


H. DONG AND H. ZHANG 


where the convergence is in the norm topology. 

We first prove an extension lemma, which is useful in our proofs. 

Lemma 2.1. Let 7 S (0,1), n>0, and f G \ (J3r^(0,ri) U (0, —r 2 ))). 

Then there exists a function F G such that 

^lR2\(Bri(0.ri)UBr2(0,-r2)) = / 

and 

< c\\f\\ n,7;R2\(Brj ( 0 ,ri)UBr 2 ( 0 .-r 2 )) > 

where C is independent of f. 


Proof. It suffices to consider the extension in 5^(0,ri). From [51 Theorem 2.19] 
and [ll Theorem 9.3], for any n > 0 there exists a solution / to the equation 

(-A)"+i/ = 0 ini?,,(0,ri), 

f = f,Dlf = Dlf,...,D:f = D:f on clS.,( 0 ,ri), 


where v is the unit normal vector of clSn (0, ri). Moreover, 

n 

ll/lln,7;Bri(0,ri) A C —i,7;aBr^ (0,7-1) ^ C'll/lln,7;R^\(Bri (0,'ri)UBr2 (0, —7-2)) J 

1=0 

where C is independent of /. Similarly, the extension of / to —r 2 ) is denoted 

by /. Finally we define 

F = /XB,.i(0,j-i) + fXBr^{0-r2) + /XR2\(Bri(0,r-i)uB,.2(0,-r2))- 

It is easy to see that F is the desired function. □ 


Denote 


_ oo - 1 , , 

a — —— e (-1,1). 

flo + 1 


Let a{x) be defined as in (O. As mentioned in the introduction, Li and Vogelius 
|9] constructed a sequence solutions to 


Di{a{x)Diu) = 0 


( 2 . 1 ) 


in R^, whose linear span is dense in iJgymda; : jxj = i?o}) for sufficiently large i?o- 
Following [^, we dehne d'j as follows: 




■E' 


ao + 1 ^ {kz + iy 


in {z : |z — i| < 1 }, 


^jiz) = (-l)" = 'i2^ +^ 0 '= 

fc=i 


{kz + iy {kz — iy 


in {2 : |z + i| > 1 and |2 — i| > 1 }, 




■E' 


ao + 1 ^ {kz- iy 


in {z : |z + i| < 1 } 
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for j odd, and 




ao + 1 ^' ' (fcz + i)^ 


in {z : I 2 ; — i| < 1 }, 


^,{z) = {-lY/^z^ +Y,{-af 


k=l 


_{kz + i)i {kz — iY 


in {z : |z + i| > 1 and |z — i| > 1 }, 


9 ^ yj 

^ (z) =-^ 

ao + 1 ^ {kz-iY 


in {z : |z + i| < 1 


for j even. Let Uj{xi,X 2 ) = i?Q-^Re'I'j(z). It is shown in [51 Proposition 8.2] 
that {uj} are solutions to (12.11) . In the lemma below, we first give an explicit 
representation of u^’s on {|a;| = Rq} in terms of trigonometric polynomials for 
Rq > 2. 


Lemma 2.2. Let Rq > 2 be a constant. For any j > 0, we have 


U2j + i{xi,X2)\\x\=Ro 


= (- 1 )^ sin( 2 j + l) 0 - 2 ^^ 
fc=lz=o 


a'^{^yy){-iy sin( 2 ; + 1)0 


( 2 . 2 ) 


for any j > 1, we have 


U2j {xi , X 2 ) I |x|=Ro 


= (-l)Jcos 2 j 0 + 2 EE 


(-a)'=(Wi)(-l)'cos( 210 ) 


(2.3) 


and 


1 — a 1 


Uo = 


1 + a oo 

Proof. From the definition, for any j >0 and x G I 80 , we have 
■U2j+i(xi,a;2) = i?(C^^'^’^^^Re4'2j+i(z) 


= Ry^+^^Re \ i-iY+Hz^^+^ + y, 


^ 2^+1 


, 2 i+l 


fc=l 


L(fcz + i)27 + l (fcz-i)27 + l 


U2j{xi,X2) = i?o ^'^Re 4 ' 2 j(z) 


= i?-^^Re (-l)^z 2 ^+^(-a)^ 

I k=l 


. 2 i 


. 2 j 


L(fcz + i) 2 i (fcz-i) 2 fJr 


(2.4) 
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Set z = i?oe'^ and we have 

U2j+l{xi, X2 )\\x\=Ro 

{ OO 

fe=i 


pi(2i+i)e 


pi(2i+i)e 


L (fci?oe'® + (fci?oe'® - i)2J+i 

1 N 2i+i 


= (-lFsi„(2j + l)« + ReX:a‘(*.ifor®«)[( ) _( ) 


k^l 


2i+l 


(2.5) 


where Ak = i/(fci?oe'^). It is obvious that jAfej < 1 for any fc > 1, and 




I + 2j 


1 N 2 i+l 
~ Ak' 


E 

;=o 


I + 2j 

2j 




fc ’ 


where f™) = .. —rr for m> n. Therefore, 

\n/ n\{m—n)] — ’ 


\\-\-Ak' Vl —Afe/ 


2 j + l 




1=0 


Plugging the formula above into (12.51) . we get (12.21) for j > 0. In the same way, for 
j > I, we obtain (1231) . Finally, (12.41) follows from a simple calculation. The lemma 
is proved. □ 


Set 

e 2 j-i = (-1)-’“^ sin(2j - 1)6», e 2 j = (-!)■’ cos2j6» 
for j > I, and Cq = uq- Clearly, for any s > 0, {ej}^g forms an orthogonal basis of 
H^y^{{x : \x\ = i?o}). It is easily seen that for any s > 0, Uj G i?syjjj({|a;| = i?o}) 
and 

ll'Wj||H«({|a;|=fl(,}) ^ + s)'*. 

Define the Hilbert space 

OO 

F = (a = (ao,ai, ■ • ■) : E 1^(1 + 

J=o 


1/2 

with the norm ||a||i« = + j)^^) ■ Then up to a constant factor, 

^sym({l^l = -^o}) is isometric to F. Denote / G F to be the infinite column vector 
ifj)j>o- Let M be an infinite dimensional matrix such that its jth column is Uj. 

Since ito = eo, we have Mq.q = 1 and = 0 for j > 1. By Lemma [121 
M = id + B, where id is the identity matrix, and B is defined as follows: for 

^J> I 


OO 

B2l-l,2j-l = -2^ 

k=l 


{kRoy+^^-^ ’ 


B2l,2j 


OO 


2E 


{kRoy+^^ ’ 


( 2 . 6 ) 




0,2j 


^■ = 2E 


(-«)^ 


B2l-l,2j — B2l,2j-1 — 0, Bio — 0. 
The following observation is crucial in our proof. 


Lemma 2.3. When Rq > 2, the infinite dimensional matrix is diago¬ 

nally dominant by column. 
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Proof. Since M = id + B, it suffices to show that J — 1- We 

first consider odd number columns. When a G (0,1), obviously B 2 i-i, 2 j-i < 0. On 
the other hand, when a € (—1,0), by (12.6p we have B 2 i-i, 2 j-i > 0. Thus, for j > 1 
we have 


GO oo fc/2Z+2j'—3 


\B21-1,2j-i \ = XI 2 XI 


a 


(2l-\-2j-3\ 
\ 27-2 J 


1^1 

Note that 


/=1 k 


oo fc/2Z+2j—3\ oo 

2W ^ ^ ) _2y-_ 


OO k f2l-\-2j — 3 


= 2XX 


a 


l2l+2j-3\ 
V 2j-2 } 


^ {Rok)^^+^l-^ I ^ ^ (i?ofc)2'+2j-2 


OO /2Z+2j —3'\ 

_^ V 2j-2 ) 


21-1 


= y _^_y 


fc=l ' 1=0 

~ fc /, I 


2j-2 + l 
23-2 


V kRf) / V kRa / 


= 11 


iklbfl- 


4(t 


)' -(t 

5:»i(yTr-‘-(uy) 


1 


fc=i 


-l/(fci?o)- 
2i-i 


l/(fci?o) 


2j-l 


2i-i 


Since lal < 1 and Ro > 2, we have 


fe=i 


1 


2i-i 


fci?0 + 1 

2i-i 


<i;((5y)""-(yT) 

k—l ^ 


= 1 . 


Then combining (12.7^ - 112.9L we obtain 

OO 

i^2/-i,2i-ii < 1. 


Z=1 


Similarly, for j > 1 we compute 


Xl-®2i.2il < XX^ 


OO OO I ^ I ^2Z+2_7 — 


Q 


Z=1 


Z=1 fc=l 




oo I It- OO /'2Z+2 j —1 

\ar v^g l 2J--1 


2Z 




oo /' 2 Z+ 2_7 —1\ 

2i-l ) 


\l=0 


=x 


a 


{kRoY 

1 


- 1 


(TTTyT)W( 


1 


2j 




oo y 

Eh‘(( 


fc=i 


^ A:i?o — 1 ^ 


1 + l/(fci?o) 

2j . / 1 \ 2i 

kRo 


- 2 


(2.7) 


( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 
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Note that for any fc > 1, i?o > 2, and j > 1, by convexity, 


Ui?o - 1/ Ui?o + 1/ ^ ■ 


Therefore, the right-hand side of (12.111) is less than 


-1 O 


1 


fc=i 


kRo -f 1 


-2f-^ 

\kRn 


which is decreasing with respect to Rq because the right-hand side of ()2.10p is 
decreasing. Thus, the left-hand side of (12.lip with |a| = 1 is less than 






Therefore, 


1^2/,2j I < 1- 

z=i 


The lemma is proved. □ 

In [9l Proposition 8.5], it is proved that for Rq sufficiently large depending on 
s, span{uj|{2;:|2;|=fl:o}} is dense in i/g®yjjj({|a;| = i?o}). In the following proposition, 
by using Lemma [2T3l we prove that i?o > 2 is sufficient to show that {uj\[x-.\x\= Ro}} 
forms a Schauder basis in = Rq}) for any s > 0. 

Proposition 2.4. For any Rq > 2 and s > 0, {Mj|{x:|x|=iJo}} forms a Schauder 
basis in H^yj^{{\x\ = i?o}); *-e-; for any f G = i?o})) there exists a 

unique sequence {aj}JLg in R such that f = OjUj. Moreover, we have 

OO 

<C(s,i?o)||/fH.({|.|=K„}). (2.12) 

i=0 

Before proving Proposition 12.41 we first show that the matrix M defines a 
bounded and invertible operator on I® in the following lemmas. 

Lemma 2.5. The matrix M defines a hounded operator on P. 

Proof. For any nonnegative integer N, we first estimate the block 
Clearly, 


C» OO OO 

X! X! "" X! y] i^2/,2ii + y] y] i52/-i,2i-ii- (2.13) 

1>N j^O 21>N j^O 21-1>N j^l 
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From ()2.6I1 and the fact that Bi^ = 0 for Z > 0, we have 


21>N j=0 21>N j=l 

°° °° \^,\kf2l+2j-l 


E 




^EE^E 


a 


rr^) 


21>N j=l fe=l 
oo 


< 


2l>Nk=l 

f 


1 \2I+1 


EEw‘ (is^) -( 


1 \2I+1 


fci?o — 1 


kRo + 1 


D 

21>N ^0 k=l 


1 


x2i+l / i \ 

) ~ \k+l/Rn) 


k — 1 /i?o 


2i+l 


fc + 1 / i?o ' 


Since Rq > 2, 


£i-i-((F^)“-(Kkn 


< 


k 1 / Rq j 




where C{a) only depends on a. Then we get 


E ElSaal < E £ CWR,-"-' 


21>N i=0 

Similarly by (1^ . we obtain 

OO 

Y E l^2Z-l,2i-l| 

21-1>N j=l 


p2Z+l 
21>N 0 


2i-l>Afc=l \ u u I / 

Combining (I2.13|) . (12.141) . and (12.151) . we have 


1>N j^O 


-N-1 

0 


(2.14) 


(2.15) 


(2.16) 


Now we are ready to show that M is a bounded operator on For any / £ t® 
with 




Z=0 


we have 


WMfwi = Yi^ + = E(i + E 

1^0 1^0 j^O 

Since Mij = Sij + by the Cauchy-Schwarz inequality 

oo 2 2 ^2 

I E I =\fi+Y^^^^M ^2 (f? + I ^ I ) . 


(2.17) 


(2.18) 


3=0 


3=0 


3=0 
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From ()2.17|1 and (j2.18[) . we get 


1=0 j=0 

Notice that (12.1611 implies for ^ > 0, 


J2lBul<CR, 

i-0 

Combining with fact that |1/||^^ < 1, we obtain 
^ 2 


-l-l 

0 


(2.19) 


^(1 + 0^*1 ^ E(i + o^*(E 


1=0 j=0 

oo 


1=0 


1=0 


+ <C{a,s), 


1=0 


□ 


where C(a, s) depends on a and s. Therefore, the proof is completed. 

Lemma 2.6. The operator on P defined by M is invertible. 

Proof. Let > 1 be a large integer to be chosen later. First we estimate the block 
of where I € [1, N], j > N: 

N [N/2] + 1 [Ar/2] + l 

EEi-^'ji- E Ei-®2i,2ii+ E E \B2i-i,2j-i\- 

1=1 j>N 1=1 2j>N 1=1 2j-l>N 

Using (12.111) . the first summation on the right-hand side of the inequality above is 
bounded by 


EE 


2|a 


k [N/2] + 1 f2l+2j-l\ 


E 


21-1 


(kRo)^i ^ (kRuY^ 

2j>Nk=l ^ 1=1 ' 


2|a| 


oo ^2Z+2j —1^ 


E 


^ ^ ^ (kRoV^ (kRo) 


2j>N k=l 


- 1 


1 


kRo -I- 1 


)"F2( 1 

/ V kRo J 


< CRo 


-N-l 


where C = C{a). Similarly, by (12.7|) and (12.8|) we have 

[A/2] + 1 

E E \B2i-i,2j-i\<CRf^-\ 

1^1 2j-l>N 

Therefore, combining (|2.20p and f|2.2ip we obtain 

E E ^ cRf^-\ 

1=1 j>N 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 


Next, we consider the N dimensional matrix Q := Since by Lemma 


2.31 is diagonally dominant by column, Q is diagonally dominant by 
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column as well, which implies that Q is invertible. We estimate Q ^ as follows. 
For odd number columns 2Z — 1 S [1, A^], by (12.811 . 


\Q2l-1.2l-l\ ~ ^ \Q2j-l,2l-l\ > \Q2l-l,2l-l\ — ^ \^hj-l,2l-l\ 

{t : 2j-lG[l,Ar] 

1 \2i-l / 1 


> 1 - 


where C{a,Ro) is a constant depending on a and Rq but not on N. For even 
number columns, by (12.111) we obtain the same estimate, 

\Q2l,2l\— ^ \Q2j,2l\ > C{a, Rq) > 0. 

{j-.2jGl2,N],j^l} 

Therefore, by [H Corollary 1] we get 

\\Q-%,<C{a,RQ)-\ 

which implies 

IIQ-^llu <C(a,i?o)■W(l + iV)^ (2.23) 

where ||Q~^||[^ is the operator norm of Q~^ in a finite dimensional subspace of P. 
Indeed, for any x = {xi,X 2 ^ • ■ - ^ we compute 

2 / \ 2 




N / N 




N / N 


\\Q-^^\\l = E E (1+< (1+E E 

1=1 \j=i j 1=1 \j=i 

Note that 

N / N \'^ / N N f ^ 

E E‘3zT/^i - E I E‘3zTi^J ^ \\Q~^\\h E 


(2.24) 


By the Cauchy-Schwarz inequality, we have 

2 
N 


i=i 


N 


i=i 


,i=i 


(2.25) 


<N\\Q-^l\\x\\l. 

Therefore, combining (I2.24p - (l2.25p . (I2.23I) is proved. 

Thanks to the estimates of the blocks, we are ready to prove the invertibility of 
M. For any hxed y G F, we need to find a unique a; € F such that Mx = y. Similar 
to F we dehne the space 

OO 

F = {a = (ai,a 2 ,---) : E < oo} 


i=i 

with the norm ||a||p = + j)^*) • Notice that the summation index 

runs from 1 instead of 0 as in the definition of 1®. 

By setting x = {xo,x)’^ and y = (?/o,y)^i where Xo,yo G K, Mx = y is written 
as 


1 (^o,i)i>i 

0 iMij)ij>i 


Xq 

X 


yo 

y 
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Clearly, it is sufficient to prove that for any y G I®, there exists x G I® such that 
j>ix = y. In fact, if there exists x G I® so that = y, it remains to 

set Xq = yo—{Boj)j>iX. From (|2.19l) . Xq is well defined. Obviously, = y 

can be further written as 

( ^ \ f \ \ 

\ B 2 id + B J \ X2 J \y2 ) ' 

where x = (xi,X2)^,y = (j/i,y2)^, xi^yi are N dimensional vectors, 

B\ = B2 — j')i>A',je[l,A] 7 B = {Bi j')i jy^. 

We rewrite the equation above as follows 

Qxi + i?iX2 = yi, i?2a;i + a;2 + .6x2 = y2- 
We take Q~^ on the hrst equation to get 

Xi = -Q~^BiX2 + Q~^yi: X2 = -B2X1 - Bx2 + 2 / 2 - 

It is sufficient to find a unique fixed point for the map on the right-hand side of the 
system above. We claim when N is sufficiently large the map is a contraction, i.e., 
the operator norm of 

/ 0 -Q-^Bi \ 

-B 2 -B J 

is less than 1. Indeed, fix x = (xi,X2)^ G I® with ||x||p = I, and we compute 

\\Q~^BiX2\\l + \\B2X1 + 5 X 211 ? 

N / N 00 \ ^ 

= E(i + ^)'^ EQfcj( E +E(i+^)'' 

k^l \j=l J k>N 

By the definition of the norm, we have 

N / N 00 \ 00 2 

^ 5,,x,) +5^(l + fc)2®(5^5,,,x,) 

k^l \i^l J k>N j^l 

( N ^ 2 \ 2 

^(l + A:)2®( ^ B,yxi) +^(l-GA:f®(^|i?,j) 

k^l / k>N j^l 

N 00 00 2 

< + E +E(i + ^)''(Ei^fc.^i) ■ ( 2 . 26 ) 

fe=l l=N+l k>N j=l 

By (12.221) . (12.191) . and the fact that \Bk,i\ < 1 for any fc, Z > 1, we know that 

N 00 N CO 

E E ^m^E E 

fc=l/=iV+l fc=l/=iV+l 



Y.\Bk,,\<CR^’^-\ 

i=i 


and 
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Therefore, combining the two inequalities above with ()2.26|1 and (I2.23|) . we have 
\\Q-^B^X2\\1 + \\B2X,+BX2\\1 
< CN{1 + + C ^ (1 + < 1, 

k>N 

provided that N is sufficiently large only depending on a, i?o, and s. Hence M and 
thus M are one-to-one and onto, and M is an invertible operator. Hence, we finish 
the proof. □ 

Proof of Proposition 12-41 From Lemmas 12.51 and 12.61 M induces a bounded and 
invertible map T on = Rq}). For any / G iLgyjjj({|a;| = i?o}), let 

T~^{f) = g and suppose that g = Since T is bounded, we have 

oo oo 

f = T[Y.^je,)=Y,a,T{e^). (2.27) 

i=o j=o 

Therefore, {uj = T{ej)} is a Schauder basis in = i?o})- By the Parseval 

identity, 

OO 

= C'(s,-Ro)||ff||ff«({|a;|=Ro}) - l^(®>'^o)||/||ff«({|x|=Ro})’ 
i=o 

which gives (12.121) . The proposition is proved. □ 


Remark 2.7. In the previous proposition, we only consider functions even in xi- 
For functions odd in xi, the same result can be proved and we only provide a sketch 


here. Dehne as follows: 


^jiz) = 




Oo + i ^ ' ' {kz + iy 


in {|z- i| < 1}, 


^,{z) = {-iy^+^y^iz^+J2{-ay 


k=l 


. (kz + i)J [kz — i)^ - 


{z : |z -I- i| > 1 and |2 — i| > 1}, 


^ (z) =- -Yi-a) 

^ ' ao -f 1 ^ 


fc =0 


{kz — i):* 


in {|2-hi| < 1} 


for j odd, and 


^j{^) = 


E- 


Oo + 1 {kz + i)f 

OO 


in {Iz- i| < 1}, 


fc=i 


{kz + iy {kz — i)P 


{z:\z + i\>l and |z — i| > 1}, 


§ (z) =- y 


ao + 1 ^ {kz- iy 


in {|z-hi| < 1} 
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for j even. Let v{x) = j{z) and following the same argument in Proposi¬ 

tion we can show that {vj} solves the equation (j4.15l) and forms a Schauder 
basis of functions odd in Xi on the circle {\x\ = Rq} provided that Rq > 2. 

We finish this section by stating [9l Proposition 8.4]. 

Proposition 2.8. Given any multi-index m there exists a constant Cm, indepen¬ 
dent of j and Rq, so that the functions Uj satisfy 

\D"^Ujixi,X2)\ < CmRo^U + 

in each of the three regions Bi fi IBi, Bi fi 052, and Bi fi 05o. 


3. Proof of Theorem o 
Now we are ready to prove our first main theorem. 


Proof of Theorem \1.1[ As explained in Remark 12.71 we only need to consider g 
even in xi. Without loss of generality, we may assume that g is smooth on dBfig = 
{|a;| = Rq}. If not, we may simply choose 2 < Rq < Rq such that K C Bn'^. By 
the elliptic regularity, M|{|a;|=Rj,} is smooth. Therefore, we can replace Rq by Rq. 
From Proposition 12.41 we have 

OO OO 

9 = '^9jUj, '^\gj\^{j + <oo (3.1) 

j=o j=0 


for some s > 3/2. It is easily seen that 



< C(j + 1 ) 


in each of the three subdomains {z : jz -|- i| < 1}, {z : |z — i| < 1}, and {z : jzj < 
i?o,|-z + i| > IJ'Z — i| > 1}, with constant C depending on Rq. Therefore, 


IWjWn^iBRg) < C{j -\- 1). (3.2) 

We consider Uk = J2 ^=q 9j'^j- Let hk be the trace of Uk on {jxj = i?o}, he., 
hk = Uk\{\x\=Ro}- Ly Proposition 12.41 in particular (12.2711 . we have hk ^ g in 
HyH{\x\=RQ}). 

Moreover, by dSH) and (EH), we get 

CO 

sup||t/fe||//i(B„^) < Cj2\9j\{j + ^) 

^ 3=0 

OO 1/2 1/2 

<c(j^lapu + if) + <„. 

3=0 3=0 

From the construction of Uj , Uk is the solution of the following equation 

Dj{aDjUk) = 0 inRi^o, Uk = hk on dB^g. 

Since hk —>■ 5 in H^^^{{\x\ = i?o}), we have that Uk ^ u in H^l^Bug). Furthermore, 
for any multi-index m, by the interior elliptic estimates, we have the pointwise 
convergence 


D^Uk{x) ^ D'^u{x) 


( 3 . 3 ) 
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with X € Bag but not on {\x — (0, 1)| = 1} and {|x — (0, — 1)| = 1}. By Proposition 
12.81 we get 

\D"^Uj{xi,X 2 )\ < CmRo\j + |to|)I™I 

in each of the three regions: Bi fi 58i, Bi fi 552, and Bi fi 55o. Therefore, by the 
Cauchy-Schwarz inequality, 

k k 

\D"^U,{x)\ < ^ \9,\\D^u,{x)\ < ^ \9j\Ro\j + \Mr 

j=0 j=0 

oo -j ^2 1 /Q 

+ (3.4) 

7=0 7=0 

for each multi-index m in each of the three regions above. From (13.31) and (13.41) . 
it follows immediately u has the desired smoothness in {|a;| < 1}. In particular, 
D'^u{x) has the same limit at the origin, whether we approach through the left 
cusp or through the right cusp. For x £ K but outside {|a:| < 1}, the piecewise 
smoothness of u follows from the classical elliptic regularity results; see, for instance, 
see [ini Proposition 1.4]. The theorem is proved. □ 

4. Non-homogeneous equations with non-symmetric coefficients 

In this section, we consider non-homogeneous equations with non-symmetric 
coefficients 

Di{a{x)Diu) = Difi, (4.1) 

where a{x) is equal to oq in i?ri(0, ri), bo in Br^{Q, —r^), and 1 in \ {Br^{Q, ri) U 

43 ^ 2 ( 0 ,—T 2 )), and ao,&o > 0. The proof is divided into three steps. We shall first 
consider homogeneous equations with ri = r 2 = 1 in Section 14.11 and then non- 
homogeneous equations with ri = r 2 = 1 in Section 14.21 and finally the general 
case in Section IT31 

4.1. Homogeneous equations. In this case, we basically adapt the proofs in Li 
and Vogelius [9], where they considered the special case oq = &o > 0. 

Recall that we use *81 and *82 to denote i3i(0,l) and i3i(0, —1), respectively, 
and *80 := R^\ (55i U $ 2 ). Let V be an open bounded subset of R^. The conformal 
mapping 2 —ijz maps 55i to {Re 2 ; > |}, *82 to (Re 2 ; < — 5 }, and *80 to (Rez G 
(—i, i)}. This leads us to study the following homogeneous equation: 



Lu := Di{A{x)Diu) = 0, 

(4.2) 

where 

A{x) = aoX{a,i>i} + + boX{xi<-^} 


Let 

Oo — 1 bo — 1 

^ 1 1 ’ ^ h 1 1 ' 

(20 “1“ 1 Oq + 1 


Choose 

a holomorphic function ^ : C \ (0, 0) —C satisfying 



II 

(4.3) 

and 

<(’(- 2 ) < when Rez > i. 

(4.4) 
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where 0 < 7 < |a/3| ^ and C is a constant. We define $ as follows: 


00 

$(z) = (1 — a) (j){z + 2k) — P(j){—z — {2k + 1)) in |Rez > 2 }’ 

k^O 

00 

$( 2 ) = (l){z) + ^ ^{aP)^[(l){z + 2k) + (j){z — 2k)) 

k^l 

- {a(l){-z + 2k - 1) + P(j){-z - {2k - 1))) in |Re 2 ; G (-i, i)|, 

OO 

$(z) = (1 — /3) ^(a/3)^’ (j){z — 2k) — a4>{—z + 2fc + 1) in |Rez < 


k=0 


Similar to [HI Proposition 8.2], we have the following proposition. 

Proposition 4.1. The function u{xi^X 2 ) = i?e$(i/z) satisfies 

Di{a{x)Diu) = 0 in K.^. (4-5) 


Moreover, u is even in xi. 

Proof. The symmetry of u in xi follows from (|4.3|) . By the property of the conformal 
mapping z —>• i/z, it suffices to verify that Re$(z) satisfies (14.211 . It is obvious that 
Re 4 >( 2 :) is harmonic in each of the three strips {xi < ^}, {xi G (—^, 5 )}, and 
{xi > 4}. It remains to check the compatibility condition. Namely, Re$(z) and 
A(a:)I?iRe$(z) are continues across the lines {xi = and {xi = —^}. Because 
4>(z) is holomorphic, by the Cauchy-Riemann equation, it suffices to verify that 
Re$(z) and A{x)D21x0. ^{z) are continues, which is equivalent to the continuities 
of Re4>(z) and A(a;)Im$(z). We only present the calculation associated with the 
continuities across the line {xi = 5 }. The verification for the case Xi = —^ follows 
the same. 

On one hand, we first compute 


00 

= (j){z) -I- (j){z) + ^ j^(a/3)^ {(fiz + 2k) + (j){z + 2k) + 4>{z — 2k) 

fe=i 

-I- (j){z — 2k)') — {a(3)^~^ {a(fi{—z -|- 2A: — I) -|- a(j){—z + 2k — 1) 
+ j3(j){—z — 2k + 1) + [3(j){—z — 2k + 1)) . 

Since z -I- z = I, the right-hand side of the equality is equal to 


(I — a) ^(q;/3)^ ((j){z + 2k) + (f){z + 2k)) 


k=0 


- P{(l){-z - {2k + 1)) + <f{-z - {2k + I))) 


which is exactly equal to 2Re<h(z)|^^^ 1 +. Therefore, the continuity of Re<I’(z) 
across the fine {xi = 4} is proved. It remains to check that ^(x)Im$ is continuous 
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and we do so by calculating 

OO 

= <l^{z) - (t>iz) + E[( a/3)^ + 2A:) — (j){z + 2k) + (j){z — 2k) — 4>{z — 2k)) 

k=l 

— {a(3)^~^ (a(^(— z + 2fc — 1) — a(j){—z + 2k — 1) 

+ l34>{—z — 2/c + 1) — P(l){—z — 2k + 1)) 

OO 

= (1 + a) ^(a/3)^ (j)[z + 2k) — (j){z + 2k) 
k=0 

— I3(j){—z — {2k + 1 )) + P4>{—z — {2k + 1)) 

because ao(l — a) = (1 + a). This completes the proof of the proposition. □ 

When 4>j{z) = l/z-i for j > 0, which is holomorphic in C \ (0,0) and satisfies 
(14.41) . from Proposition l4.ll Uj := i?()’^Re4'j(z) is a solution to (14.51) for each j with 


= (1 — a) y^(a/3)^' Y. - ——^ + Py -—- in IBi, 

^ ^ L(i + 2A:z)J '^(i+(2A: + l)z)jJ 

OO 


fe=i 


- (i + 2kzy (i — 2kzy 


+ /3i 


(-i + (2fc + l)z)l (i + (2A: + l)z)l. 


in So, 


4-,(z) = (l-/3)^(a/3)^ 

for j odd, and 


fe=0 


.(i —2fcz)f (—i + (2fc + l)z)f- 


= (1 - a)^{al3Y 


k=0 


. (i + 2fcz)l (i + {2k + 1 )^)-^ 


in S2 


in Si, 


J J 


— a 


-Ptz- 


{—i+{2k + l)zy (i + (2A: + l)z)J- 


in So 


4i,(z) = (l-/3)^(a/3)^ 


— a 


k=0 


-(i —2fcz)f (—i + (2fc + l)z)f- 


in S2 


for j even. By straightforward calculations, we have the following propositions 
similar to [U Propositions 8.3 and 8.4]. 

Proposition 4.2. Given any integer m, there exists a constant Cm, independent 
of j, so that the functions ipj satisfy 
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in each of the three regions 

[z : |z| < l,|z + i| > l,|z-i| > l|, 

: |z| < 1, |z — i| < l|, and : |z| < 1, |z + i| < l|. 

Proposition 4.3. Given any multi-index m, there exists a constant Cm indepen¬ 
dent of j and Rq , such that the functions Uj satisfy 

\D"^U,{x,,X 2)\ < CmRo\j + hi)'™' 

in each of the three regions: Bi fl Bi fl and Bi fl ®o- 


Next, we investigate Uj restricted on {|x| = Rq} with Rq > 2. By setting 
z = for j > 0, 


U 2 J +1 = (-1)-^ sin(2j + 1)0 


Re 


CX) 

{E(»«‘ 


oi(2i+l)e 


pi( 2 i+l)e 


fc=l 


L(i + 2A:i?oe*®)2J+i (i - 2fcRoe*®)2J+i 


A( 2 i+i)e 


A( 2 i+i)e 


(-i + (2/c + (i + (2fc + l)i?oe'®)2f+i 

U 2 j = (-!)■’ cos(2j0) + Re I y)(a/3)‘ 


,i 2 ie 


fe=i 


L(i + 2fci?oe‘®)2i (i-2fcRoe*®)2J 


ni 2 ie 


— a 




.i 2 te 


(-i + {2k + l)i?oe‘®)2f (i + {2k + l)Roe^O)^i 


r]}- 


It is easy to see that 

(a-l)(/3-l) 

° 1-0/3 ’ 

U 2 j = (-!)■’ cos(2j6>) + 0{Rq'^^) on {|a;| = i?o}, j > 1, 

M 2 j+i = (-1)^ sin(2j + 1)6» + on {|x| = i?o}, 1 > 0. 

Therefore, similar to [HI Propositions 8.5 and 8.6], we obtain the following denseness 
result on {uj}. 


Proposition 4.4. Given any s > 0, there exists a constant Cg < oo, so that 
span{uj\i\x\=Ro}} *5 dense in Hly^{{\x\ = i?o}) provided that Rq > Cg- Moreover, 
for any function g G Hgymi{\x\ = .Ro}); we may approximate it by g = liuikhk, 
with 

Nk 

hk =^"/j''''uA{\x\=Ro} e 
j=o 


and 

k 1 /Q 

j=0 


where C depends on s and Rq, but independent of k and g, and Nk G N. 
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Remark 4.5. In Proposition 14.41 the solutions are even in xi. For the case when 
solutions are odd in xi, we define $ as follows: 


l>(z) = (1 - a) ^(a/3)'= U{z + 2k) + - {2k + 1)) 




{Rez > -}; 


$(z) = 4>{z) + ^2 + 2fc) + <i){z — 2k)) 


+ (a/3)'=-i {a(t){-z + 2fc + 1) + Pc^{-z - {2k + 1))) 


in {Rez £ (- 


1 1 
2’2 


)}; 


$(z) = {1 — P) ^^(a/3)^ (j){z — 2k) + a(j){—z + 2k + 1) 


k=0 


{Rez < 


and v{xi,X 2 ) = Im4)(i/z). It is easy to see that v is odd in xi and following the 
proof of Proposition l4.ll we can show that u is a solution to (14.51) . Further, we can 
prove similar result to Proposition 14.41 


Now we are in a good position to show the following theorem, the proof of which 
is similar to that of Theorem 11.11 and follows the lines in proving 0 Proposition 
8 . 1 ]. 

Theorem 4.6. Suppose Rq > Co, where Co is a constant depending on ao and bo- 
Let g be in H^^‘^{^{\x\ = i?o}); O'lT'd u € H^{Bjig) denote the weak solution to 

Di{a{x)Diu) = 0 in Bfig, u = g on {jxj = i?o}- 

Then 

u G C°°(R:\ (*Bi U*B 2 )), uGC°°($i), and uGC'°“($ 2 ) 

for any compact set K C Rrq . 


Proof. Without loss of generality, we assume that g is smooth and even in xi as 
in the proof of Theorem 11.11 Let h^ be the approximating sequence of g as in 
Proposition 14.41 with some fixed s > 3/2, i.e., 


Nk 

hk^g in i4®({|a;| = i?o}) and hk 

3=0 


with 

/ \ 1 /2 
(El7fV(l+Jr) <C\\g\\H^, 

3=0 

where C depends on Rq. By straightforward calculations, we have 



< C{J + 1 ) 


in each of the three subdomains {|z + i| < 1}, {|z: — i| < 1}, and {jzj < Ro, |z + i| > 
1,1-2^ — i| > 1}, where C depends on Rq. Hence, 


IWjWmiBRg) < C{j + 1). 
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By the Cauchy-Schwarz inequality, the sums Uk = convergent in 

H^{Bro) and 


Nk 

j=0 

/ \ 1 /2 / \ 1 /2 
<c(E 17^1^0-+ 1)^0 + <-■ 

j=0 j=0 

By the linearity of the equation, Uk is the solution to 

D,{a{x)D,Uk) = 0 inBijo, Uk = hk on {|a;| = i?o}- 

From our construction, hk ^ g in = Rq})- Thus we have Uk ^ u in 

H^{Br^). From the elliptic regularity theory, we know that for any multi-index m, 

D'^Uk{x) D'^u{x) (4.6) 

at any point inside Sflo, but not on the circles {|x—(0, ±1)1 = 1}. From Proposition 
oi we get that for any multi-index m 


Nk Nk 

|i7™C/fe(a;)| <^|7f ||i7™«,(x)| IVO’ + IH)'™' 

j^O j^O 

/ \ 1 /2 / \ 1 /2 

(4-7) 

j=o j=o 

in each of the three regions i?i fl fBi, Si fl I 82 , and Si fl IBo- From (14.61) and 
it follows immediately that u has the desired smoothness properties in Si. 
For X G K but not in Si, the piecewise smoothness of u follows from the classical 
elliptic regularity results. □ 


4.2. Non-homogeneous equations. In this subsection, we consider the non- 
homogeneous equations by constructing Green’s function of the operator in (14.11) . 
By applying the conformal mapping z —>■ i/z, we shall first construct Green’s func¬ 
tion of the operator L defined in (14.21) . i.e., 

Di{A{x)DiG{x, y)) = S{x - y), 
where Di is with respect to Xi. 

Denote k = (fc, 0), where fc G Z. Let y = {yi, — 2 / 2 )- It is well know that 
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In this section, for simplicity of exposition, we write A log \x — y\ = 5{x — y). We 

.1 1 ] 

2 ’ 2 ^ 


define G{x,y) as follows: when yi € {—h, h), 


G{x,y) = {l-a)'^{al3f log |a; + 2k - ?/| 

-^log|a; + 2k + l + y| in{xi>i}, 

OO 

G{x, y) = log lx - 2/1 + ^ [(a/3)'' (log |x + 2k - i/l + log |x - 2k - y\) 




- (a/3)'' ^ (/3 log |x + 2k - 1 + 2/1 + a log |x - (2k - 1) + 2 /|) 
in {xi e (-i, i)}, 

OO 

G(x, 2 /) = (1-/3)^(a/3)'' log |x-2k- 2/1 


- alog|x - (2k + 1) + 2/1 in{xi<--}; 


when 2/1 > 


G{x, 2 /) = log |x - 2/1 + a log |x - 1 + y\ 

~ log |a; + 2k + 1 + 2 /| in {xi > i}. 




G(x, 2 /) = (l + a)^(a/3)'' log |x-2k-//I 




1 1 . 


-/31og|x + 2k+1 + 2/1 in {xi e (--,-)}, 

log k - 2k - 2/1 in {xi < -i}; 




when 2/1 < — 5 , 


G{x,y) = log |x + 2k - y\ in {xi > i}, 


1 4- ao 




G(x, 2 /) = (l + ^)^(a/3)'' log|x + 2 k- 2 /| 




1 1 . 


-alog|x- (2k + l) + 2 /|J in {xi e (--,-)}, 

G(x, y) = log |x - 2/1 + /3 log |x + 1 + y\ 

- log |x - (2k+l) + y| in {xi < . 

° fc =0 

Proposition 4.7. The function G{x,y) defined above is Green’s function of L, i.e., 

Di{A{x)DiG{x, y)) = S{x - y). 
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Proof. To show G{x,y) is Green’s function, it is sufficient to prove that for j/ € 
a.e., AG(x, y) = S{x — y) for x ^ {xi = i}U{a;i = and G{x, y), A{x)DiG{x, y) 
are continuous in x across the two lines {xi = ±^}. 

We first consider the case when 2 /i G (— 5 , 5 )- It is obvious that for xi < —^ 
and k > 0 


|a;-k-?/|> 0 , |a;-k + ?/|> 0 , 

which implies AxG{x, y) = 0. Similarly, we can check that for xi > ^, AxG{x, y) = 
0. Moreover, for xi G (—1/2,1/2), note that 

Aj, log|x -yl = S(x - y), 


and 


Aj; log \x + 2k-y\= Aj, log |x - 2 k - y\ 

= Aj, log \x+ 2k -l + y\ = A^ log |x - ( 2 k - 1 ) + y| = 0 


provided that fc > 1. Thus, 


AxG{x,y) 


S{x-y) forxe(-i,i). 


It remains to verify the continuities of G{x,y) and A{x)DiG{x,y) across the lines 
{xi = 1/2}. For simplicity, we only present the calculations associated with the 
case xi = 1/2. We first check that G{x,y) is continuous at Xi = 1/2. By a 
straightforward calculation, we have 


G(x,y)|^^^i- = 2log 


. 00 

2 


- yif + (x2 - yif 


1 


/c=l 


log {{2k + - - yif + (x 2 - y 2 f) 


+ log {{2k - - + yif + (x 2 - y 2 f) 

- {aPf~^ /31og {{2k - 1 + i + y{f + (x 2 - y 2 f) 
+ a log (( 2 fc - 1 - i + (x 2 - y 2 f) } 

^00 . 

’ log {{2k +-- yif + (x 2 - y 2 f') 

fc =0 


- /3 log {{2.k +'^+yif + {X 2 - 2 / 2 )^) 
= G{x,y)\^^^, + . 
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Next we check that A{x)DiG{x,y) is continuous across {xi = 1/2}. We compute 
A{x)DiG{x,y)\^^^i~ = Di\og\x - yW^^^i 

CO 

+ Di log |2k + a; - i + £>i log |2k - a; + 


k^l 


- PDi log |2k - 1 + a; + - aDi log |2k - 1 - a; - 

2 fc + 1 — j/i 2 fc + I + ?/i 


L(2fc+1 - yi)2 + (a;2 - 2 / 2 )^ (2A: + | + yi)^ + (x2 - ^ 2 )^ 

On the other hand, we have 

A{x)DiG{x,y)\^^^i+ = aoDiG{x,y)\^^^i 

CO 

= ao(l - a) y^(a/3)^ 


2fc + 1 - 


/c =0 


L(2fc+l-yi)2 + (a;2-2/2)^ 


-/3 


2 /c + I + 2/1 


Since 
we get 


{2k + f + 2 / 1 )^ + {x2 - yiY 
ao(l - a) = (1 + a) 


A{x)DiG{x,y)\^^^i- = A{x)DiG{x,y)\^^^i + . 

When yi > 1/2, the singularity appears in the region {xi > 1/2}. For complete¬ 
ness, we present the calculations below. We first verify 

AxG{x,y) = 5{x - y). 

For xi > 1/2 and k > —1, it is easy to see that |k -|- a; -|-y| >0, which implies 

Aa, log |a; -f k + y| = 0 . 

Combining with the fact that 

A2,log|x - 2/1 = 5{x - 2/), 


we get 

A^G{x,y) = 6{x - y) for xi > 1/2. 

Similarly, for xi G (—1/2,1/2), 

|x — 2/1 > 0, |x-I-2k — 1-I-2/1 > 0, and |x —2k —2/|>0 
provided that fc > 1. Thus, 

A^G(x,y) = 0 for xi G (-1/2,1/2). 

In the same way, we have 

Aj;G{x,y) = 0 for xi < —1/2. 

Next we verify the continuities of G{x,y) and A{x)DxG{x,y) at xi = ±1/2. By 
the same argument as in the case 2/1 G (—1/2,1/2), without loss of generality, we 
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only check the continuities at xi = 1/2. To this end, we compute 

1 OO - 

y)Li=i- = —^ log [{2k - - + yif + {X 2 - y 2 f) 


/c =0 


l3log [{2k + ^ + yi)2 + {x 2 - y 2 f 


1 + a 


log ((2 - Vif + {^2 - y 2 ?) 


+ (1 + a){a/3 - I3)Y^ log [{2k + ^ + yi)^ + {x 2 - y 2 f 


k=0 


On the other side of xi = 1/2, we calculate 
1 + a . 


G{x,y)\^i+ = 


Since 


■iog((^ -yif + (^2 -y^f) 


^'*'1 + “) f: 1^ log ((1 + m + 2 *. + 1 )^ + (X 2 - »)/ 


1 + Oo 


fc =0 


{1 + a){al3 - 13) = - 


2/3(1 + ct) 


1 + Oo 

it follows immediately that 

G{x,y)\^^^i~ = G{x,y)\^^^, + . 

Next, we verify that A{x)DiG{x,y) is continuous at Xi = 1/2 and compute 

OO 

A{x)DiG{x,y)\^^^i- = (l + a)y](a/3)'= Hi log |2k - a; + 

- /3Hi log |2k + 1 + a; 

2k-\+yi 


= -{^ + a)'^{al 3 ) 

-/3 


k^O 


l-( 2 fc- i +?/l )2 + (X 2 - ?/ 2 )^ 


2 fc + I + ?/i 


( 2 A:+|+?/i )2 + (x 2 -y 2 )^ 


= -(1 + a ) 


yi - 2 


(yi - + ix 2 - y 2 )‘^ 


-(l + Q;)2/3^(a/3)'‘ 


2 fc + yi + I 


/c =0 


(2fc+f+yi)2 + (a;2-?/2)^ 


On the other hand, we have 

A{x)DiG{x,y)\^^^i+ = -ao(l - a!)yY 


yi- 2 


{h-yi)^ + {x 2 -y 2 )^ 


2ao/3(l + a) 

1 + Oo 


Em)V 


h 3~ yi 2k + 1 


k—O 


{k+yi + 2k + 1)2 + (a ;2 - 7 / 2 )^ 


Because 


ao(l — a) = (1 + a) and (1 + a)^/3 = 


2ao/3(l + a) 

1 + Oo 























ON AN ELLIPTIC EQUATION ARISING FROM COMPOSITE MATERIALS 


27 


it follows that 


A{x)DiG[x,y)\^^^i- 


A{x)DiG{x,y)\^^^^+. 


Therefore, G{x, y) satisfies 

D^{A{x)D^G{x, y)) = 5(a; - y) (4.8) 

for the case yi > 1/2. Similarly, we can check that (14.8p holds when j/i < —1/2 as 
well. The details are omitted. Thus, G{x, y) is Green’s function to the divergence 
type operator L and we complete the proof of the proposition. □ 


Now, let us turn back to the original operator Lu = Di(a{x)Diu{x)). The 
conformal mapping z ^ ijz can be written in real variables as 0 : —?► 


0 i(a;) 


X2 

2 I 2 ’ 

Xi + X2 


02 (x) 


Xl 

^1 + ^2 ' 


For any integer fc, denote Xj^[x) = 0(0(a:) + k), which is a conformal map. 
According to G, we define G{x,y) as follows: when y G IBo, 


G{x,y) = (1 - [log\X 2 kix) - y\ - /31og |A 2 fe+i(a:) - y\] for x G *Bi, 

k =0 

OO 

G{x, y) = log \x-y\G^ ^(a/3)'’’ log |A 2 fe(a;) -y\+ log \X_ 2 k{.x) - 

k^l 

- p\og\X 2 k-i{x) -y\ - alog|A_( 2 fc_i)(x) -y\ for x G *Bo, 

OO 

G{x,y) = (l-/3)^(a/3)'’’[log|A_2fc(a:)-2/| 

- a\og\X_(^2k+i){x)-y\] fora;ie*B2; 

when ?/ G *Bi, 

G{x, y) = log lx - ?/| + a log I A_i (x) - y\ 

_ 2^(1^^ ^ fc for xG*Bi, 

1 + On 


fc=0 


G{x,y) = (l + Q:)^(a/3)'’’[log|A_2fc(x) - y\ - p\og\X 2 k+i{x) - y\] for x G *Bo, 

k—0 

G{x, y) = y'(a/3)'' log |A_ 2 fc(x) - y\ for xG‘=B 2 ] 

1 + oo 


/c=0 
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when y G * 82 , 


G{x, y) = log \X 2 k{x) - y\ for xi G Si, 

1 + ao ^ 


A;=0 
00 


G(x,y) = (1 + log \X 2 k{x) - y\ 

- alog|X_( 2 fc+i)(x)-y| for X G *Bq, 
G(x,y) =log|x-y|+/31og|Xi(x)-y| 

2a(l+/3) 


1 + 


■ log |X_2fe-i(x) - y| forxG®2. 


fe =0 


Proposition 4.8. The function G(x,y) defined above is Green’s function of L. 


Proof. Similar to the verification of G(x, y) being Green’s function of L, we first 
check that 


Aj;G{x, y) = S(x — y) for x ^ 9(iBi U * 82 )- 

In order to show this, we consider the case when y G $0 as an example. When 
X G $1 U * 82 , we show that G(x, y) is harmonic. For instance, when x G i 8 i, 

0(x) + 2k G {xi > 1 / 2 } and 0(y) G {xi G (—1/2,1/2)}. 

Therefore, 0(x) + 2k 7 ^ 0(y), implying X 2 k{x) 7 ^ y. Similarly, we have X 2 k-i{x) 7 ^ 
y. Combining with the facts that A^, log |x — y| =0 when x 7 ^ y in and that Xk 
is conformal, we obtain that G{x,y) is harmonic in $ 1 . In the same way, we can 
show that G{x,y) is harmonic in $2 as well. When x G $ 0 , as we mentioned in 
the beginning of this section, we use the notation log |x — y| = 5{x — y). Each 
term in the expression of G{x,y), with the exception of log |x — y|, is harmonic in 
*80 by the same argument in proving AxG{x,y) = 0 in * 81 . Hence, when x G I 80 , 
Aa;G(x, y) = 5{x — y). For the case when y G *81 U * 82 , the same argument can be 
implemented to show that AxG{x,y) = 5(x — y), and we omit the details. 

It remains to verify the continuities of G(x, y) and a(x)D^G{x, y) across the two 
circles {|x — (0, ±1) = 1|}, where v is the unit normal vector field of cl(i 8 i U * 82 ). 

Because 0 is a conformal map, the continuities of G(x,y) and a(x)D^G{x,y) 
is equivalent to the continuities of G(0(x),y) and a(0(x))ZIi(G(0(x),y)), respec¬ 
tively. Note that a(0(x)) = A{x) and 

^ \y\ 

| 0 (y)| | 0 (x)|’ 

which by similarity of triangles implies that 

|Q(a^) -1/1 ^ M 

|x- 0 (y)| |xr 


(4.9) 
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We take the case when y S as an example and the other cases can be verified 
in the same way. Since 0(0(a;)) = x, we have 

G(0(a;),y) = log|0(ai) - y\ +alog|0(a;- 1) -y\ 

- log |Q(a^ + 2k + 1) - y| in {xi > i}, 

1 + ao ^ 




^(0(2;),^) = (1 + a) log |0(2: - 2k) - y\ 


/c=0 


1 1. 


-/31og|0(a; + 2 k+l)-y|J in {xi e 
G(0(a;),y) = ^(a/3)'= log |0(x - 2k) - y| in {a;i < -i}. 




By taking (03 into account, G{Q{x),y) has the expression: in {ii > ^} 


G(0(a:),y) = log |a: - 0(y)| +alog|x- 1 + 0(y)| 

— y^(Q:/3)^ log |a; + 2k + 1 + 0(y)| — log |x| — a log \x — 1| 


1 + Oo 

2^(1 + q) 

1 + Oq 


fc=0 

oo 


log |a; + 2k + 1| + log l^li 


fe=0 


1 — aP 


in {si e (-5,5)} 

OO 

G(0(a:),y) = (1 + a) ^^(a/?)'' log |a; - 2k - 0(y)| - ^log |a; + 2k + 1 + 0(y)| 




- (1 + a) y](a/3)'' log |a: - 2k| - /31og |a; + 2k + 1| 


fe =0 


(l + a)(l-/3) 
1 — a/3 


log|y|; 


and in {xi < —i} 

G(0(a;),y) = log\x - 2k - e{y)\ 

1 + On ^ 




2(1+ OL, , (l-/3)(l + a) 


1 + &o 


■y](a/3)'"log|2:-2k| 


fc=0 


1 — aP 


log|2/|. 


Observe that 


G(0(x),y) = G(x,0(y)) - H{x) + log |y|, 


where G{x, y) is Green’s function of L and H{x) is a function obtained by replacing 
y with 0 in the expression of G{x, y). Since we verify the continuities of G{x, y) and 
A{x)DiG{x, y) across the lines {cci = ±i} in Proposition l4.71 the same proof shows 
that H{x) and A{x)DiH{x) are continuous across the lines {cci = ±i}. Therefore, 
G{0{x),y) and A{x)DiG{0{x),y) are continuous by linearity, and G{x,y) is the 
desired Green’s function. □ 
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With the help of Green’s function constructed above, we are ready to consider 
the non-homogeneous equation 

Di{a{x)Diu{x)) = Difi{x) 

in general 2? C Now, we state our theorem in the case when ri = r2 = 1. 

Theorem 4.9. Let n € N U {0} and 7 £ (0,1). Assume that u is a weak solution 
to the equation 

Di{a{x)Diu{x)) = Difi in V, 

where 

a(x) = ao mlBi, a(x) = bo in ^ 82 , a(x) = 1 m *80, 
and for each i, fi is pieeewise C'"”''', i.e., 

Let Vg = {x € V, dist(x, dV) > e} for e > 0. Then u is pieeewise in Vg up 

to the boundary, i.e., 

u £ n »i), u £ n <82), u £ n »o)- 

Proof. We prove the theorem by considering two cases. 

Case 1: (0,0) ^ V. Define Lli = V D 18i,r22 =1)0 I82, and LIq = V O Sq- 
Because (0,0) ^ V, it is easy to see that any point in 1) belongs to at most two 
subdomains Lli, which is exactly the case in [3 Remark 3(ii)]. Therefore, we apply 
[II Theorem 2 and Remark 3(ii)] to obtain that when n = 0, u £ piecewise in 
2?e, i.e., for any 0 < z < 2, 

uGC^''^{Venni). 

For n > 0, we use an induction argument. For a ball away from the circles 
{|a; — (0,±1)| = 1}, the conclusion follows from the classical Schauder estimate 
for Poisson’s equation. We only need to consider a ball Bi{x) C Ve and x £ 
{|a; — (0,±1)| = 1}. Notice that by locally flattening the boundary, it is sufficient 
to consider 

DiioijDjU^ — 

where are piecewise smooth in and R?. with bounded derivatives and hi £ 
(7ii,7(R^), g C"'’'^(R^), where R^ (or R?_) is the set of points on the plane such 
that X 2 > 0 (or X 2 < 0). Here we only give a sketch of the proof. For n = 1, by 
taking derivative with respect to the tangential variable xi, we have 

Bii^aij Dj Diu') — DiDihi Dii^D^ipXij'^Djuf 

Thanks to the case n = 0, we have that Du £ C'*' piecewise, which implies the right- 
hand side can be written as Difi for some piecewise functions fi. Therefore, we 
apply the result of the case fc = 0 to obtain that DDiu £ C'*' piecewise. It remains 
to estimate D^u, which is obtained from the formula 

= af^iDihi — DifuijDju) — D2{a2iDiu)). 

Therefore, DD 2 U is piecewise as well. By induction, it is easy to prove that 
u £ piecewise. 

Case 2: (0, 0) £ V. There exists 0 < / < 1 such that Bi C V. We define a cutoff 
function rj £ C^{Bi), which equals 1 on R//2. Let v = up, which satisfies 

DifaDivfx)) = Diffip + auDip) - fiD^p + aD^uDip 


(4.10) 
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in . We define 

u{x) = - [ Dy.G{x, y)fi{y) + G{x, y){fi{y)Di'q{y) - aDiuDiy) dy, (4.11) 
where 

fi = fiV + auDif]. (4-12) 

Since G is Green’s function of L, the function u defined above is a solution to (14.101) . 
When u is restricted to {|x| > |}, the result follows from Case 1. Therefore, it 
remains to estimate u\bi/ 2 - Since v = u in B 1 / 2 , it suffices to consider v instead of 
u. Because 

Di{a{x)Di{u - u)) = 0 

in by Theorem 14.61 with a sufficiently large we know that u — v is piecewise 

smooth. Hence, it suffices to estimate u instead of v. 

Note that supp(Ili?7) C {1/2 < \x\ < /}, on which by Case 1 it G (^"+1,7 
piecewise. Combining with the definition of fi in (I4.12p . we get that fi are piecewise 
C"’^. By (IdUll . 

u{x) = - Dy.G{x,y)fi{y)dy - G{x,y){f^{y)D,r]{y) - aD^uDiy) dy 

:= ui{x) + U 2 {x). 

Since the estimates of ui and U 2 are quite similar, we only consider ui as follows: 


ui{x)=- Dy^G{x,y)f,{y)dy 


/R2 


= - / Dy.G{x,y)f^{y)dy - / Dy.G{x,y)fi{y) dy - / Dy,G{x,y)f,{y) dy 

J<S2 J<So 

:= -wi{x) - W 2 {x) - W 3 {x). 


We focus on the case a: G Ho H i?i and the same argument can be applied to the 
other cases as well. By the definition of G{x,y), we have 


where 


wi{x) = (1 + a) V(a/?) 




Dyi log \X_ 2 k{x) - y\My) dy 


Si 



Dy, log \X 2 k+i{x) - y\My) dy'j 


= {l + a) [h{X. 2 k{x)) - ph{X 2 k+i{x ))], 


(4.13) 


h{x)= Dy,\og\x-y\f,{y)dy, h{x)= Dy, log \x - y\My) dy. (4.14) 




'<81 


Since log \x — y\ is the fundamental solution of the Laplace equation in h(x) 
satisfies 

Ah =-Di{fiX<Bi) inR^. 

Since / is piecewise C'”’''" and the interface 9*81 is smooth by the same method in 
dealing with Case 1, we obtain that h G C'"“''^’'’'(ili) piecewise. 








32 


H. DONG AND H. ZHANG 


For any /c 0, by the definition, 


= Re 


, Im ■ 


= ^ ( Re 


1 


kz + i 


7 ,1 + Im 


1 


kz + i 


i/z + k' i/z+ 

where z = Xi + \X 2 - Notice that for a; G * 80 , jfcz + i| > 1. By a straightforward 
calculation, it is easily seen that for any j > 0 and a; G i?i fl So, 


\D^Xk{x)\<C\k\^+\ (4.15) 

where C is independent of k. Moreover, since x G flonBi, Q{x) G {|a;| > l}n{a;i G 
(— 5 , 5 )}, which implies that 0(a;) — 2k G {a: : Xi < 1/2, \x\ > 1}, and 

X_ 2 k(x) = 0(0(x) - 2k)) G Ri \ ^. (4.16) 

Therefore, combining (14.151) and (I4.16L with the chain rule, we have for any k > 0 


\\h{X_2k{xmn+i,r,B,nno < C{k + 1)^+^ 



It remains to estimate h{X 2 k-i)- Notice that 


/25i 


-Dy,log|a;-^/|/^(?/)dy = / Hy, log |a; 


l<S 2 


implying that 

Ah{x) = D,{f,{x)x<S 2 )- 
Note that for any k > 1 and a; G Ri D IBq, 


y\My)dy , 


X 2 k-i{x) = 0(0(a;) + 2k - 1) G Ri n *Bi. 

Furthermore, the regularity of fi{y) is the same as fi{y), which indicates that 
h{X 2 k-i) can be estimated in a similar way as h{X- 2 k)- Hence, combining the 
estimate of h and h, we have 


00 

< C(1 + a) y^(a/3)^(||fe(X_2fc)|| n+l,7;Bin*Bo + ^ || ^(^2fc+l) || n+l,7;Bin*Bo ) 
k^O 
00 

l)"'''^(ll^l|n+l,7;SinB2 + ll^l|n+l,7;Bin<Bo + II ^l|n+l, 7 ;Bin'Bi ) 

k—0 

2 

- ^^2 (ll/i^lU.7;®i + ll^^i’llUoi'Bi) + ^{\\Hl^{B2) + ll^llioo(B2)) 

i=i 

2 

- ^(^^2 ll/ilU.7;®l + ll“lln,7;®iU(Bi\B,/2))i (4-17) 

where the last term on the right-hand side is estimated in Case 1. Here in the last 
inequality above, we use the fact that 

II^IIAoo(B2) + II^IILoo(B2) 

2 2 

< c'22 ll/ilUoo('Bi) < + II^I|Loo(B,\B,/2))> 

i=i i=i 

which can be deduced directly from (14.141) . 
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By symmetry, it is easy to see that the argument applied to wi can be imple¬ 
mented to W 2 as well. We omit the details. Now let us estimate W 3 with a little 
modification. Similar to the expression in (I4.13L we have 

00 

Wsix) = g{x) + ^(a/ 3 )'' {g{X2k{x)) + g{X_2k[x)) 

k=l 

+ Pg{X 2 k+i{x)) + ag{X_^^ 2 k+i){x))^ , 

where 

g{x)= DyAog\x-y\fi{y)dy, g{x) = - Dy Jog \x - y\f^{y) dy. 

J<Bo J<So 

Since fi{y) has the same regularity as fi, it is sufficient to consider g, which satisfies 

Ag = -A(/zX<Bo)- 

Here, we cannot directly apply the result in Case 1 because of the singularity of 
the domain Sq. Nonetheless, by Lemma [ 2.11 there exists Fi € which is 

the extension of fiX<Bo and satisfies 

ll^i||n.7;R2 < C (||/i || + l|w||n,7;(BASi/2)n<Bo)' 

Therefore, define 


9 


/R2 


Dyi log \X 


y\Fjy)dy, 


01 := 



Dyi log|a: - y\Fjy) dy, 


and 


02 := 


/»2 


Dy. log \x 


y\Fi{y)dy, 


which satisfy 

Ag = -DiFi, Afli = -DjFiX^sJ, and Aga = -A(-F*X' 82 )- 

From the classical Schauder estimate, we have g € By the estimate of 

wi above, we have gi, g 2 are all piecewise (7”+^’'*', which implies 5 = g — 0 i — g 2 is 
piecewise as well. Then we can follow the same argument in the estimate 

of wi (cf. (14.171) 1 and obtain a similar estimate for 1 x 3 

2 

1^311 n+l,7;Bin^0 — c(x:ii/jik 7;BinBo + lklln,7;(Bi\Bi/2)nSo)- 
Hence, we show that u is piecewise and the proof is completed. □ 


4.3. Two balls with different radii. Next, we consider the general case that the 
two balls have different radii. Specifically, 

a{x) = aoXB,.j( 0 ,ri) + ^OXB,.^ ( 0 ,-r- 2 ) + Xx>\(B,.^(o,rl)uBr 2 (o,-r• 2 ))■ 

Denote dBr^ (0, ri) = Cl and (0, —r 2 ) = C 2 . By scaling and reflection, without 
loss of generality, we may assume that r 2 > ri > 1/2. Now, we consider a conformal 
map T : C —>■ C,T{z) = where zq € C and zq = zi -|- Z 2 i. It is well known 

that for zo ^ Cl U C 2 , T maps Ci and C 2 to two circles. We shall find a suitable zq 
such that the two circles have the same radius. Indeed, by a simple computation it 
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is easy to see that T{Ci) and T{C 2 ) are circles with radii \ri/{zl + — 2riZ2)\ and 

\’i' 2 l(z\ + + 2 r2Z2)\, respectively. Then, we only need to find (zi, Z 2 ) such that 

|ri/(zi + z| - 2riZ2)\ = \r2/{zf + z| + 2r2Z2)\. 

It is obvious that 


ri{z'^ + z| + 2r2Z2) = ?'2(zi + z| - 2riZ2) 

(?'2 - ri)izl + zl) = 4ziZ2?'ir2. (4.18) 


Note that (I4.18p can be written as 

£1 _l_ _ 4 rir 2 

Z2 Zl r2 - ri 


(4.19) 


which implies the existence of (zi,Z 2 )- Moreover, if (zi,Z 2 ) is a solution of (I4.19L 
for any s > 0 (szi,sz 2 ) solves (14.191) as well. Hence, we can pick zq outside V, 
which indicates that T is smooth in V. After a translation, rotation, and scaling 
of the coordinates, we can assume that T maps i?ri(0,7’i) and 5 ^ 2 ( 0 ,— 7 - 2 ) to *81 
and * 82 , respectively. Therefore, we obtain the desired conformal map T, which is 
a diffeomorphism between T) and T{'D). 

Now we are ready to prove Theorem 11.11 


Proof of Theorem 11.11 In order to consider the equation 

Di{a{x)Diu{x)) = Difi{x), 

we define v(x) = u(T~^x) in TifD). By a straightforward calculation, we have 
DMT~\x))D,v) = Dk{mh{T-^x)) 

in T{'D), where r]ki are some smooth functions in T[T>). The operator on the 
left-hand side of the equation above is the same as the operator in Theorem 14.91 
Therefore, applying Theorem 14.91 we get v S piecewise in T{'D) C\Bi, where 

Bi is a small ball around the origin. This, combined with the smoothness of T and 
the chain rule, implies that u G piecewise in a small ball B^ C T~^{Bi) 

around the origin. For the region outside B-^, it follows from Case 1 in the proof of 
Theorem 14.91 Hence, the theorem is proved. □ 


By a standard freezing coefficients argument, finally we prove Corollary 11.31 

Proof of Corollarv \1.3[ First, we consider a ball Bi{xo) such that Bi{xq) C Pe D 
-Sri(0,ri). By the classical Schauder estimate, we have 

2 

ll'*^l|n-|-l,7;Bi/2(a:o) — ^ II'*^11 Aoo(Bi (xq)) T ^ ^ II/j II n,7;Bi (xp) V (4.20) 

i=i 

Similarly, we can show (14.201) holds if Bi{xq) C (~) 5 ^ 2 ( 0 ,—r 2 ) or Bi C Vg \ 
{Br^{0 ,ri)U Br^{0,-r2)). 

Second, we consider a ball i?;(xo) such that xq G 9 ( 5^1 (0, ri) U Hj. 2 (0,—r 2 )) but 
0 ^ Bi{xq). By flattening the boundary, this is essentially the same as Case 1 in 
the proof of Theorem 14.91 This implies that u is piecewise in this case. 
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It remains to consider Bi{xo) = 5/(0,0) C Vg. In this case, we define a as 
follows: 

a{x) = lim a{x), x € Brj(0,ri), 

x^{0,Q) 
x^Br^ (0,ri) 

a(x) = lim a(x), x € Br^iO,—r 2 ), 

a;->(0,0) 
x^Br-2 (0,— 1 - 2 ) 

a(x) = lim a(x), otherwise. 

a:->(0,0) 

a;GR^\(Brj (0,1-1)05^2 (0, —r2)) 

Then the equation can be written as 

Di{aDiu) = Di{fi + (a — a)Diu). (4-21) 

Denote fli = Bi ni?ri(0,ri), CI 2 = Bi fl (0, —r 2 ), and CIq = Bi\ {Brj^{0,ri) U 

-Sr 2 (0i’’ 2 ))- Let = Ofc n Bi /2 for 0 < fc < 2. By the Leibniz rule and the 
interpolation inequality, we get 

[(a - + {P + 

where e > 0 and C(e) depends on e,n, and ||a||„_.y;n,^. 

Applying Theorem 1 1.21 to (I4.21L we obtain 
2 

'y ||w||n+l, 7 ;f 2 fc 
k=0 

2 2 2 2 

^ ^(")( E E + E II“IIl.(o.)) + C,{P + £) EM„+i.,;a,, (4.22) 

k—Qj—1 fc—0 fc—0 

where Ci is independent of e. Since 

Mn+l,7;nfc — C'Mn+l.Tififc + L*)!!] , 

where C is independent of u and I, by taking I and e sufficiently small in (I4.22L we 
obtain 

2 22 2 2 

E \W\\n+l,T,^k - ^(E E ll■('tlln,7;nfc + E ll“llLoo(nfc)) + ^ E [“]n+1.7;S^fc\f2fc ’ 

fc—0 fc—0 i— 1 fc—0 fc—0 

where is estimated in the previous case. Therefore, thanks to the 

argument of partition of the unity, the corollary is proved. □ 
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